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Abstract: Super symmetry breaking in a metastable vacuum is re-examined in a cosmological context. 
It is shown that thermal effects generically drive the Universe to the metastable minimum even if it 
begins in the super symmetry-preserving one. This is a generic feature of the ISS models of metastable 
supersymmetry breaking due to the fact that SUSY preserving vacua contain fewer light degrees of 
freedom than the metastable ground state at the origin. These models of metastable SUSY breaking 
are thus placed on an equal footing with the more usual dynamical SUSY breaking scenarios. 



1. Introduction 



It has recently been noted that dynamical supersymmetry breaking in a long-lived metastable vacuum 
is an interesting alternative to the usual dynamical supersymmetry breaking scenario [l] 1 . It was shown 
that a microscopic asymptotically free theory can, for certain choices of parameters, be described by 
a Seiberg dual theory in the infra-red which has a metastable minimum at the origin. Moreover, the 
global supersymmetric minima are also located where the macroscopic theory is well under control, 
and it could be shown that the tunnelling rate to the true global minima could be made parametrically 
small. This class of theories has recently attracted a lot of interest [3-14]. It opens up new avenues 
for SUSY phenomenology, but one open question is that of naturalness; can one explain the fact that 
the Universe today resides in a metastable minimum? Are these models more, or less, natural than 
the usual scenario of dynamical supersymmetry breaking? 

In this paper we point out that ending up in such a metastable minimum is in fact very natural and 
even generic in the context of a thermal Universe. The theory at the metastable minimum has (quite 
generally) more massless and nearly massless degrees of freedom than at the supersymmetric minimum, 
an unusual consequence of the fact that the SUSY preserving vacua appear non-perturbatively 2 in the 
Intriligator-Sciberg-Shih (ISS) model. We will show that due to these relatively light states, at high 
temperatures (greater than the supersymmetry breaking scale), the metastable minimum at the origin 
becomes the global minimum. Moreover, we will see that even if the fields begin in the supersymmetric 
minimum, high enough temperatures will thermally drive them to the SUSY breaking minimum where 
they remain trapped as the Universe cools. More precisely, we consider a scenario where, at the end 
of inflation, the Universe is very cold and we may assume that it is in the energetically preferred 
supersymmetric vacuum. Then the Universe reheats. If the temperature is high enough it will then 
automatically evolve towards the SUSY breaking state as we will show in the following. Thus these 
models offer an alternative and appealing explanation for why supersymmetry is broken: the early 
Universe was driven to a supersymmetry breaking, metastable minimum by thermal effects. 3 

In general we think of the settings where the ISS model forms a sector of the full theory which 
includes a supersymmetric Standard Model (MSSM). Supersymmetry of the full theory is broken if 
the ISS sector ends up in the metastable non-supersymmetric vacuum. For the full theory to be driven 
to this SUSY breaking vacuum by thermal effects we, of course, have to assume that the relevant fields 
of the supersymmetry breaking sector (ISS), of the MSSM sector, and of the messenger sector were at 
some time in thermal equilibrium. This is the case if the SUSY breaking scenario is gauge mediation, 
direct mediation, or even a visible sector breaking. (On the other hand, a SUSY breaking sector which 
couples to the MSSM only gravitationally would remain out of thermal equilibrium.) 

'in a simple modified O'Raifeartaigh model the possibility of SUSY breaking in a metastable vacuum has been 
considered already a while ago in [2]. It is interesting to note that one of the motivations for this work was that the 
Universe may get stuck in such a metastable state. 

2 In other words, supersymmetry is restored dynamically. 

3 If the initial conditions are such that the Universe starts in the supersymmetry breaking state (the case considered 
in Refs. [15,16] subsequent to this paper), it will stay there (cf. also [1,17]). 
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The fact that the Universe can be driven to a metastable minimum by thermal effects is well 
known in the context of charge and colour breaking minima in the MSSM [18, 19]. However these 
models differ from our case in two respects. First, it is in the global minimum away from the origin 
that a symmetry, namely supersymmetry, is being restored in our scenario. In addition the effective 
potential is extremely flat because the SUSY preserving true vacua are generated dynamically. It is 
therefore automatically very sensitive to finite temperature effects, whilst having a metastable vacuum 
that is extremely long lived. 

In this paper it will be sufficient to consider the ISS metastable SUSY breaking sector on its own. 
The thermal transition to the non-supersymmetric vacuum will be driven by the above mismatch 
between the numbers of light degrees of freedom in the two vacua of the ISS model. 

2. The ISS model at zero temperature 
2.1 Set-up of the model 

The Intriligator-Seiberg-Shih model [1] is described by a supersymmetric SU(N) gauge theory coupled 
to Nf flavours of chiral superfields ipf and (p % c transforming in the fundamental and the anti- fundamental 
representations of the gauge group; c = 1, . . . , N and i = 1, . . . , Nf. There is also an Nf x Nf chiral 
superfield which is a gauge singlet. The number of flavours is taken to be large, Nf > 3N, such 
that the /3-function for the gauge coupling is positive, 



the theory is free in the IR and strongly coupled in the UV where it develops a Landau pole at the 
energy-scale A^. The Wilsonian gauge coupling is 



The condition (2.1) ensures that the theory is weakly coupled at scales £ < A^, thus its low-energy 
dynamics as well as the vacuum structure is under control. In particular, this guarantees a robust 
understanding of the theory in the metastable SUSY breaking vacuum found in [1]. This is one of the 
key features of the ISS model(s). 

One notes that this formulation of the theory can only provide a low-energy effective field theory 
description due to the lack of the asymptotic freedom in (2.1), (2. 2). At energy scales of order A^ 
and above, this effective description breaks down and one should use instead a different (microscopic) 
description of the theory, assuming it exists. Fortunately - and this is the second key feature of the 
ISS construction [1] - the ultraviolet completion of this effective theory is known and is provided by 
its Seiberg dual formulation [20-22]. The microscopic description of the ISS model is N = 1 super- 
symmetric QCD with the gauge group SU (N c ) and Nf flavours of fundamental and anti-fundamental 



b = 3N - N f < 



(2.1) 




(2.2) 
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quarks Q, t and Q l . The number of colours in the microscopic theory is N c = Nf — N and the number of 
flavours Nf is the same in both descriptions. It is required to be in the range N c + l < N f < §iV c . The 
lower limit on Nf ensures that Eq. (2.1) holds, while the 5™ lcro coefficient of the microscopic theory is 
positive (the /3-function is negative and the microscopic theory is asymptotically free) 

6™ = 3N C - N f , ^N c < bf CI ° < 2N C -1 (2.3) 

This microscopic formulation of the ISS model will be referred to as the SU(N C ), Nf microscopic 
Seiberg dual. It is weakly coupled in the UV and strongly coupled in the IR. However, we already 
mentioned that the vacuum structure of the theory and the supersymmetry breaking by a metastable 
vacuum should be considered in the low-energy effective description of the theory, which is IR free. 
This description of the ISS model is known as the SU(N), Nf macroscopic Seiberg dual formulation. 
From now on we will always assume this macroscopic dual description of the ISS model. For the 
purposes of this paper, the microscopic dual description is only necessary to guarantee the existence 
of the theory above the Landau pole A^. 



The tree-level superpotential of the ISS model is given by 

W cl = hTr Nf (p^<p - V 2 Tr^$ (2.4) 

where h and fx are constants. The usual holomorphicity arguments imply that the superpotential (2.4) 
receives no corrections in perturbation theory. However, there is a non-perturbative contribution to 
the full superpotential of the theory, W = W c \ + Wd yn , which is generated dynamically. Wd yn was 
determined in [1] and is given by 

This dynamical superpotential is exact, its form is uniquely determined by the symmetries of the 
theory and it is generated by instanton-like configurations. 



The authors of [1] have studied the vacuum structure of the theory and established the existence 
of the metastable vacuum |vac) + characterised by 

ifp) = <V3 T > = A* ' <*> = ' V + = ( N f ~ (2-6) 

where V+ is the classical energy density in this vacuum. Supersymmetry is broken since > 0. In 
this vacuum the SU(N) gauge group is Higgsed by the vevs of (p and <p and the gauge degrees of 
freedom are massive with m gaug e = g\x. 

This supersymmetry breaking vacuum |vac)+ originates from the so-called rank condition, which 
implies that there are no solutions to the F-flatness equation = for the classical superpotential 
W cl in (2.4), 

F#= hititf-vPSfr^O (2.7) 
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The non-perturbative superpotential of (2.5) gives negligible contributions to the effective potential 
around this vacuum and can be ignored there. It was further argued in [1] that the vacuum (2.6) has 
no tachyonic directions, is classically stable, and quantum-mechanically is long-lived. 

In addition to the metastable SUSY breaking vacuum |vac) + , the authors of [1] have also identified 
the SUSY preserving stable vacuum 4 |vac)o, 

(<p) = (y5 T > = , <$)=$„= M7o t Nf , V = (2.8) 

where Vo is the energy density in this vacuum and 

(Nj:—3N \ -1 
he -rrj=ir\ ) and e:= iL <1 (2.9) 

This vacuum was determined in [1] by solving the F-flatness conditions for the complete superpotential 
W = W c \ + lUdyn °f the theory. In the vicinity of |vac)o the non-perturbative superpotential of (2.5) 
is essential as it alleviates the rank condition (2.7) and allows to solve F^j = equations. Thus, 

i 

the appearance of the SUSY preserving vacuum (2.8) can be interpreted in our macroscopic dual 
description as a non-perturbative or dynamical restoration of supersymmetry [1]. 

When in the next Section we put the ISS theory at high temperature it will be important to know 
the number of light and heavy degrees of freedom of the theory as we interpolate between |vac)+ and 
|vac)o- The macroscopic description of the ISS model is meaningful as long as we work at energy 
scales much smaller than the cut-off scale A^. The requirement that e = /i/A^ <C 1 is a condition 
that the gauge theory is weakly coupled at the scale fi which is the natural scale of the macroscopic 
theory. Equation (2.9) implies that 70 S> 1 or in other words that there is a natural separation of 
scales /j « $0 = ^7o *C A^ dictated by e <C 1. The masses of order h$o will be treated as heavy and 
all other masses suppressed with respect to h<&o by a positive power of e we will refer to as light. In 
the SUSY preserving vacuum |vac)o the heavy degrees of freedom are the Nf flavours of cp and <p, they 
get the tree-level masses m v = M>o via (2-4). All other degrees of freedom in |vac)o are light. 5 

In the vacuum |vac) + all of the original degrees of freedom of the ISS theory are light. Hence we 
note for future reference that the SUSY preserving vacuum |vac)o has fewer light degrees of freedom 
than the SUSY breaking vacuum |vac)+. The mismatch is given by the Nf flavours of <p and (p. The 
fact that the metastable SUSY breaking vacuum has fewer degrees of freedom than the supersymmetric 
vacuum is an important feature of the ISS model. As we will explain below, this property will give us 
a necessary condition for the thermal Universe to end up in the metastable vacuum in the first place. 

4 In fact there are precisely Nj — N = N c of such vacua differing by the phase e 2lrl /( N f~ N ) as required by the Witten 
index of the microscopic Seiberg dual formulation. 

5 Gauge degrees of freedom in |vac)o are confined, but the appropriate mass gap, given by the gaugino condensate, 
is parametrically (in e) smaller than m v = /i3>o- Thus the gauge degrees of freedom can still be counted as light for 
sufficiently high temperatures. 
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2.2 Effective potential 



Having presented the general set-up, let us now turn to the effective potential between the two vacua 
|vac) + and |vac)o of the ISS model. This can be determined in its entirety as follows. We parameterise 
the path interpolating between the two vacua (2.6) and (2.8) in field space via 



(p(a) = (p T {a) = an I N 
y UN f -N 



$( 7 ) = j/j,l Nf , 0<7<7o,l>cr>0 



(2.10) 



Since the Kahler potential in the free magnetic phase in the IR is that of the classical theory, the 
effective potential V can be determined from knowing the superpotential of the theory (2.4), (2.5). 
First we will calculate the classical potential V c \ along path (2.10) and ignoring the non-perturbative 
contribution W^ yn . Using (2.4) we find, 

^2 V d (7, °) = Tr^l^^l 2 + Tv Nf \^\ 2 + TY^Ito - /A/ = H 4 (2AT 7 V + N{a 2 - l) 2 + N f -N) 

(2.11) 

This expression for V c \(^,a) is extremized for a = y/l — 7 2 or a = 0. We substitute the first solution 
for a into for the range of < 7 < 1, and the second solution, a = for 1 < 7 < 70. We thus 
obtain 



1 



Ki(7) 



N f - N + 2iV7 2 (l - i 7 2 ) 0<7<1 
Nf 1 < 7 < 70 



(2.12) 



This potential is accurate for 7 < 1, where the potential is rising. In order to address the flat regime 
1 < 7 < 7o w e need to include the contribution of the non-perturbative Wd yn . This gives 

1 2 



V = Tr« 



N I h N f rr f — 1 -hn 2 &. 



.Nf-SN 



|ftVl^((^)- 



(2.13) 



Combining these expressions we have the final result for the effective potential interpolating between 
the two vacua: 



Nf — N + 2NY(1 - i 7 



Nf-N \ 2 

— _ 1 



< 7 < 1 



1 < 7 



(2.14) 



For our forthcoming thermal applications we have included the subscript T = to indicate that this 
effective potential is calculated at zero temperature. We plot the zero-temperature effective potential 
(2.14) in Figure 1. The key features of this effective potential are (1) the large distance between the 
two vacua, 70 3> 1, and (2) the slow rise of the potential to the left of the SUSY preserving vacuum. 
(For esthetic reasons 70 in Figure 1. is actually chosen to be rather small, 7 = 70 = 7.5.) 



We note that the matching between the two regimes at 7 = 1 in (2.14) holds up to an insignificant 

Nf-N 

l/7 n correction which is small since 70 3> 1. This tiny mismatch (which can be seen in Figure 1.) 
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Figure 1: Zero temperature effective potential Vr=o(j) of Eq. (2.14) as a function of 7 = For the SUSY 
preserving vacuum |vac)o we chose 7 = 70 = 7.5. The SUSY breaking metastable minimum |vac)+ is always at 
7 = 0, and the top of the barrier is always at 7 = 1. We have taken the minimal allowed values for ./V and Nf, 
N = 2, N f = 7. 




can be easily corrected in the derivation of V by using the full superpotential in both regimes, rather 
than neglecting Wd yn at 7 < 1. However, it will not change anything in our considerations. Other 
corrections which we have dropped from the final expression include perturbative corrections to the 
Kahler potential. These effects are also expected to be small at scales much below the cut-off Al- 

The authors of [1] have already estimated the tunnelling rate from the metastable |vac) + to the 
supersymmetric vacuum |vac)o by approximating the potential in Figure 1 in terms of a triangle. The 
action of the bounce solution in the triangular potential is of the form, 

^bounce = ^2 -ffp { — J ~ e 4(iV / -3iV)/(iV / -JV) > 1 f ° r £<1 ( 2 ' 15 ) 

Thus as argued in [1] it is always possible, by choosing sufficiently small e, to ensure that the decay 
time of the metastable vacuum to the SUSY ground state is much longer than the age of the Universe. 

On closer inspection the constraints imposed by this condition are in any case very weak. In order 
to estimate the required value of e, note that the Euclidean action ^bounce gi yes the false vacuum decay 
rate per unit volume as 

r 4 /U = D 4e - 5 b?unce. (2.16) 

D4 is the determinant coefficient and is irrelevant to the discussion. To determine whether the Universe 
could have decayed (at least once) in its lifetime, we multiply the above by the space-time volume of the 
past light-cone of the observable Universe. This results in a bound to have no decay of roughly [19,23] 

Sb?u„cc~ 400. (2.17) 



This translates into an extremely weak lower bound on 
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Note that the bound is on the relative width to height of the bump. Thus it is indeed simply the 
flatness of the potential which protects the metastable vacuum against decay. In terms of e = /jl/Al 
the bound (2.18) reads 



It is interesting to note that this very weak bound (2.18), (2.19) becomes strong when expressed in 
terms of e itself in the minimal case of N = 2 and Nf = 7, 



In the following Section we will explain why at high temperatures the Universe has ended up in the 
metastable non-supersymmetric vacuum in the first place. 



3. Dynamical evolution at finite temperature 

3.1 The shape of the effective potential at finite temperature 

The effective potential at finite temperature along the $ direction is governed by the following well- 
known expression [24]: 



The first term on the right hand side is the zero temperature value of the effective potential; for the 
ISS model we have determined it in Eq. (2.14) and in Figure 1. The second term is the purely thermal 
correction (which vanishes at T = 0) and it is determined at one-loop in perturbation theory. The 
rii denote the numbers of degrees of freedom present in the theory 6 and the summation is over all of 
these degrees of freedom. The upper sign is for bosons and the lower one for fermions. Finally, mj($) 
denote the masses of these degrees of freedom induced by the vevs of the field <I>. 

We have noted in the previous Section that as we interpolate from the metastable vacuum to the 
supersymmetric one, the Nf flavours of <p and dp acquire masses = h& = hji r ) and become heavy at 
large values of <I>. To a good approximation all other degrees of freedom can be counted as essentially 
massless. 7 As we are not interested in the overall additive (T dependent, but field independent) 
constant in the thermal potential, we need to include in (3.1) only those variables whose masses vary 

6 Weyl fermions and complex scalars each count as n — 2. 

7 In the companion paper [17] we will refine this analysis by including contributions of gauge degrees of freedom to 
(3.1). This will lead to an even lower bound for the reheat temperature, T R , necessary for the Universe to end up in the 
non-supersymmetric |vac) + . 
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(2.20) 
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Figure 2: Thermal effective potential (3.2) for different values of the temperature. Going from bottom to top, 
the red line corresponds to the temperature > C rit where we have only one vacuum at 7 = 0. The orange 
line corresponds to 6 « 6 cr it where the second vacuum appears and the classical rolling stops. The green line 
is in the interval Gdcgcn < © < @crit where one could hope to tunnel under the barrier. The blue line is at 
~ Qdogcn where the two vacua become degenerate. Finally, the black line gives the zero temperature potential 
where the non-supersymmetric vacuum at the origin becomes metastablc. 

from being light in one vacuum to being heavy in the other. These are precisely the 99's and (p's. This 
implies 8 



The above expression is written in terms of a dimensionless variable 7 = $/// and we have also defined 
a rescaled temperature 



In Figure 2 we plot the thermal potential (3.2) for a few characteristic values of temperature. 

In general one can consider a very wide range of temperatures < T -C Al corresponding to 
< B <C l/(he). From Fig. 2 where we have plotted the effective potential at various temperatures 
we immediately make out several interesting temperatures. At high enough temperature there is only 
one minimum at the origin and one would expect to be able to roll down classically to the non- 
supersymmetric minimum |vac)+, as will be discussed in Section 3.2.3. Below = cr it a second 
minimum forms. It will turn out, that for our ISS potential this critical temperature will actually be 
less than 70. In the following analysis, we will consider the range of temperatures hfi < T < h&o 
corresponding to 1 < < 70. 

At temperatures below cr it the SUSY breaking minimum at the origin still has lower free energy 
then |vac)o- We will calculate the bubble nucleation rate to the true vacuum |vac)+ in Section 3.2.2. 
Another distinguished temperature ©dcgcn is when the SUSY preserving vacuum |vac)o becomes de- 
generate with the non-supersymmetric vacuum |vac)+. For a generic potential one would expect that 

8 There are only two terms in the sum in (3.2): one for bosons (+), and one for fermions (-). The prefactor 4NNf 
counts the total number of bosonic degrees of freedom in <pf and (p % a . 




(3.2) 



e = T/\hfx\ 



(3.3) 



- 8 - 



somewhere between these two scales there is a temperature 0* (©degen < ©* < ©crit) where the rate 
of bubble nucleation turns from large to small. Our main goal is to determine the temperature where 
the transitions (classical or tunnelling ones) from |vac)o to |vac) + stop. 



3.2 Transition to a SUSY breaking vacuum and the lower bound on Tr 

In this subsection we discuss the two possibilities for the field to evolve to the SUSY breaking vacuum 
|vac) + , classical evolution and bubble nucleation. 



3.2.1 Classical evolution to the SUSY breaking vacuum: Estimate of T cr i t 



It is clear from our expression for the thermal potential that at sufficiently high temperatures there 
is only one minimum, namely |vac) + , and the SUSY preserving minimum disappears. This is caused 
by Nf flavours of (p and dp becoming heavy away from the origin. The disappearance of |vac)o is most 
easily seen analytically in the limit of very high temperatures 7^0 where the thermal potential 
grows as 2 7 2 in the 7 i.e. $ direction, 

V© - U e =o(7) ~ NN f G 2 7 2 - const 6 4 , 9 > 7 = $/// (3.4) 

Thus there is a critical temperature cr it such that at > cr it there is only one minimum, and at 
< ©crit the supersymmetry preserving second minimum starts to materialise. 

A better estimate can be obtained by using the approximation, 

± dq q 2 In (l T exp(- yY + 7 2 /© 2 )) ~ ©' (^@) § exp (-^) , for 7 » 6 (3.5) 

for the integral in the 6 dependent part of Eq. (3.2). Comparing the derivatives (first and second) of 
the 6-dependent part and the 6-independent part of the effective potential V@ in the vicinity of 70 
we get an estimate, 

70 ^_ V2ttI ( N f ^ 2 



which we have also confirmed numerically. 

Although the height of the barrier in the zero temperature potential is only of the order of fi 
where f( < 70/i = h&o the temperature necessary to erode the second minimum is only slightly 
(logarithmically) smaller than h&o due to the exponential suppression ~ exp(— 7/8) apparent in Eqs. 
(3.2) and (3.5). 

In the early Universe the situation is not static and the temperature decreases due to the expansion 
of the Universe. So one should check that even in the absence of a second minimum the field has time 
enough to evolve to |vac)+ before the temperature drops below O cr it- We will comment on this issue 
more in the next subsection. Here, we just point out that the temperature drops on a time scale 
Mpi/T 2 > 1/T. 
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3.2.2 Bubble Nucleation: Estimate of % 



Let us now turn to the possibility of bubble nucleation, with an estimate of the temperature where 
the transitions from |vac)o to |vac) + change from fast to slow. Bubble nucleation is of course only- 
possible above the temperature where the two minima become degenerate, so let us first estimate this. 
We compare the effective potential (3.2) at the origin 

Ml = 0) = (N — N f ) - (n B + 7 -n F ^j (3.7) 

to the value of the effective potential at the second vacuum. The latter remains approximately zero 
since the temperature Td ege n is much below <l>o- In (3.7) ra_B(ra_F) counts the number of bosonic(fermionic) 
degrees of freedom (ANNp in both cases). We thus find 

e ^={ i hWNN J ) y> f ° r 7o>>l (3 - 8) 

Now we would like to know when bubble nucleation is fast enough to lead to a phase transition. Thus 
we consider temperatures in the range ©degen < © < © C rit such that the second vacuum |vac)o is 
already formed but still higher than |vac)+. In the following we want to estimate 0* where tunnelling 
becomes too slow. 



For this purpose we will now derive a simple estimate on the action of the tunnelling trajectory 
in the thermal effective potential of the our model where the temperature is in the range ©degcn < 
© < ©crit- The potential we want to model is depicted in Figure 2, its characteristic features are 
that the 'false' vacuum |vac)o is far away in the field space from |vac) + and that the potential climbs 
very slowly from |vac)o to a shallow barrier and then descends steeply to the 'true' vacuum |vac)+. 
These features are a reflection of the fact that at zero temperature the vacuum |vac)o was generated 
non-perturbatively (as reviewed in the previous Section). We will model this class of potentials with 
a simple linear no-barrier potential 9 , 

Min(7) = 0{ri-i)N f K(i-Tj) (3.9) 

where K and i] are constants, Nf is put for convenience, and 6 is the step-function. We plot this 
potential in Figure 2 for a convenient choice of constants along with the exact thermal potential. It 
is clear from this figure that the tunnelling rate for our model potential will always be a little higher 
than in the real case. 



The potential (3.9) admits a simple analytic solution for the bounce configuration. The tunnelling 
in this potential was discussed in detail in Ref. [26] in the T = case. We have performed a similar 
calculation in the 3-dimensional settings relevant to the thermal case. The tunnelling configuration is 
the bounce solution 7bounce(^) which extremizes the 3-dimensional action 

±S™=±Anj drr 2 {^ Nf {<S>'f + V T (<S>)) = N f ^- 2 J dxx 2 (i( 7 ') 2 + ± V&nfr)) (3.10) 

9 We have checked the dependence on the characteristic scales of the potential using other simple approximations, for 
example two matched quadratics and a triangular potential. 
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Figure 3: A simple no-barrier model for the thermal effective potential. The modelling potential Vim (7) is 
shown in red and the exact thermal potential Vr(7)is black. The tunnelling rate to the true vacuum for this 
model is calculated in the text. 



with the appropriate initial conditions. The right hand side of (3.10) is written in terms of our usual 
dimensionless variables 7, and the rescaled radial distance x = h[ir. The classical equation for the 
bounce is then 

7 " + -y = w 9 7 %) (3.11) 

For the case (3.9) the bounce solution reads 

*iz (312) 

where x m = is the matching point for the two branches of the bounce solution and of its 

derivative. The field configuration (3.12) describes a bubble of size x m with 7 = 3r] on the outside at 
large x 3> x m , and 7 = on the inside at x = 0. The asymptotic value of 7t>ounce is the 'false' vacuum 
|vac)o, thus we should identify 2>rj with 70. 



The action (3.10) on the bounce trajectory is 

1 3D 4vr 8V6V 2 "* ,, 

f ^bounce " @p iV / -5- -JT V- L6 > 

Now using the identification 3rj = 70 and NtKr) oc AV, where AV is the drop in the potential which 
should be taken to be AV oc 4 we get 

^Sw* « const (^) 3 (3.14) 

using the original dimensional variables $0 anc i T. In order for the Universe to have undergone a 
phase transition by bubble nucleation one requires a sufficiently high nucleation probability of bubbles 
of metastable vacuum. The nucleation rate is given by 

T 3 ~ T 4 e~ s3 °/ T . (3.15) 



This rate must be integrated from a maximum (reheat) temperature Tr to the temperature ld e g e n ~ M 
at which the metastable and supersymmetric minima are degenerate. The fraction of space remaining 
in the broken phase is e~ p where [19] 



M 4 rT ' 



.,, I" T- 2 (l-T Q /Tfe~ s3D l T dT (3.16) 

•'T'dcgon 

Mi, 

^ < 3 ' 17 > 
where To is the temperature of the Universe today. This gives 

D 256-^- M W 

P ~ e R ^3D ( 3 - 18 ) 

and the bound becomes jP- < 1.3. From this we conclude that the temperature T* at which the 
tunnelling transitions are still possible is very high (at the very top of its defining range), 

T* > <*>o , or 0* > 7o (3.19) 

The temperature where bubble nucleation becomes significant has the same parametric dependence as 
the critical temperature cr it (our simple estimate is insufficient to capture logarithmic dependencies) 
and we conclude that T* ~ T cr i t . 



3.2.3 Behaviour of the field after nucleation/rolling 

If the temperature rises above the critical temperature, the field is in principle free to roll; however, is 
the critical temperature a sufficient condition for the field to always end up at the origin? Assuming 
that the Universe evolves in the standard FRW manner, we should check that the phase transition has 
time to complete, and that the falling temperature does not "overtake" the field. In order to model 
the field after the transition, we can approximate the potential at temperatures above T cr i t as linear 10 ; 

. <£ 

V = const T — const' (3.20) 

Neglecting for the moment the effect of the tp&tp coupling, the field equations are 

$ + 3H$ = -const (3.21) 

We may assume that the Hubble constant H ~ T 2 /M Pl is negligible in this equation, so, taking initial 
values T > T crit i.e. T ~ <3?o, we find that the field falls to the origin in time 

Ai roll ~ |J ~ 1 ~ M~}tl (3.22) 



Here, we want just a rough estimate of the time scale of the evolution towards the supersymmetry breaking minimum. 
For this we use a simple adiabatic approximation (i.e. we use classical evolution equations but with an effective thermal 
potential) and a simplified potential. One can argue that, in general, some kind of non-adiabatic treatment is required. 
Nevertheless, we expect that our conclusions from this simple estimate remain qualitatively correct. 
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where on the right hand side to defines the time at which the temperature is Tq, and we have considered 

1 

a radiation dominated Universe with scale factor a(t) ~ (t/to) 2 ■ The result (3.22) should be compared 
to the typical time necessary for the Universe to cool to temperatures where the second minimum 
|vac)o becomes prominent. The evolution of T is given by 

1 

Tm =W>~ T °(W (3 ' 23) 

For example the time needed to reach T = T$/\f2 is At coo i = to which gives At coo i » At ro n. This 
means that the field is able to roll essentially instantaneously to the origin where it undergoes coherent 
oscillation. 



Having concluded that the fields roll down to |vac) + sufficiently fast, we still need to check that 
the oscillations of the fields around |vac) + would not bring us back to |vac)o as the Universe cools. It 
is easy to see that if there was only the Hubble damping the $ field would oscillate out of the origin 
as the Universe cools. Indeed, the <I> oscillations essentially preserve energy in a comoving volume; 



^l^max-R 3 = const > (3.24) 

where R ~ 1 2 is the scale factor in a radiation dominated Universe, and m$ = 4^ is the temperature 
induced mass of $ at the origin. In an adiabatic regime RT = const so we find 

- & 

The field would then escape from the origin because when T ~ u, (assuming To ~ $0) the size of the 
oscillations would be 

$max ~ > fi. (3.26) 

Fortunately, any open decay channel of $ will typically provide a damping rate T$ ~ T, capturing 
the field in |vac)+. An example are couplings to the messenger sector fields, /, of the form h'f&f. 
Since we assume thermal equilibrium, such a coupling must exist and cannot be arbitrarily weak. The 
decaying oscillation amplitudes are of order 

$max(T) _ [T ' l r$(t _ to) 

Sufficient damping always occurs; imposing $ miB < fj, requires only that ~ log(<E>o/^) which is 
always satisfied and the Universe ends up in the SUSY breaking state |vac)+. Note, however, that the 
number of oscillations before damping is proportional to log(<E>o//u) and may be large. 



3.2.4 Lower bound on reheating temperature Tr 

Now, let us summarize our analysis and explain under which conditions it is natural for the early 
Universe to settle down at the metastable SUSY breaking vacuum. At the end of inflation the Universe 
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is in a state of very low temperature and one may assume that it is in the energetically preferred 
supersymmetric state 11 |vac)o- After inflation the Universe reheats to a temperature Tr. Already at 
relatively low temperatures ~ 7d eg e n ~ A* the supersymmetry breaking vacuum will have lower free 
energy than the supersymmetric vacuum. However, if Tr falls in the range ?degen < Tr < T* ~ T cr it, 
the Universe will remain in the state |vac)o although |vac) + is energetically preferred. It is stuck there 
because the barrier makes classical evolution impossible and bubble nucleation is too slow. Above 
T* the bubble nucleation rate will be high and above T cr i t the classical evolution becomes possible. 
Hence, the Universe will evolve to the preferred supersymmetry breaking state |vac) + . We conclude 
that if the reheating temperature Tr fulfills 

Tcrit ~ ^N^Nj-N) ~ Tr ~ f ( 3 - 28 ^ 

it is ensured that the Universe always ends up in the supersymmetry breaking ground state 12 |vac)+. 

In principle one should also consider the possibility of a transition back towards the SUSY vacuum 
for T < Tdegen- It is known that this does not happen at zero temperature [1] . We have made a simple 
estimate of this effect at T > and concluded that the Universe remains trapped. More recently, this 
has been discussed in depth in [15-17]. 

The condition (2.18) or (2.19) that the metastable vacuum |vac) + is long-lived at zero temper- 
ature does not necessarily require $o to be very large or € 1 - 2N /( N f- N ) to be very small. For a low 
supersymmetry breaking scale n of order of a few TeV the lower bound on the reheating temperature 
in (3.28) can be easily satisfied for reheating temperatures as low as around 10 or 100 TeV. On the 
other hand, for models with a significantly higher SUSY breaking scale, this bound becomes more 
constraining. 

4. Conclusions and discussion 

We have examined the consequences of meta-stable SUSY breaking vacua in a cosmological setting. 
As noted in [1], the ISS theory naturally admits a finite number of isolated supersymmetric vacua 
(as determined by the Witten index) along with a larger moduli space of metastable SUSY breaking 
vacua. Since the latter is a much bigger configuration space, ISS suggested that it is more favorable 
for the Universe to be populated in the metastable SUSY breaking vacua. In this paper, we have 
considered thermal effects in the ISS model and have shown that the early Universe is always driven to 
metastability: as long as the SUSY breaking sector is in thermal equilibrium this provides a generic, 
dynamical explanation why supersymmetry became broken. This phenomenon is a consequence of 
two distinguishing properties of the ISS theory which are not necessarily shared by other models 
with metastable SUSY breaking vacua, namely, that the metastable minima have more light degrees 
of freedom than the SUSY preserving vacua, and that the metastable vacua are separated from the 

n If it is already in the supersymmetry breaking vacuum, |vac) + , it will stay there. 
12 In this estimate we ignore weak logarithmic corrections in e. 



- 14 - 



SUSY preserving vacuum by a very shallow dynamically induced potential. Both play significant roles 
in our argument and our various estimates. 



Given the current interest in "landscapes" of one variety or another, it is difficult to resist specu- 
lating on their implications for the ISS models, given the conclusions of our study. One longstanding 
problem that can be addressed in this context is that of the hierarchy between the Planck and su- 
persymmetry breaking scales. Consider the possibility that the SUSY breaking sector has in fact a 
large number of product groups, with a landscape of metastable ISS minima into which the Universe 
could be thermally driven, with a range of values of SUSY breaking. The very existence of the lower 
bound on Tr implies that there is a maximal value of a SUSY breaking scale n which characterizes 
metastable vacua we can reach as the Universe cools down. Using our estimate for Tr, Eq. (3.28) we 



where Al is the Landau pole of the theory. Thus not only is the SUSY breaking determined to lie 
below Al, it is also parametrically suppressed by powers of Tr/ Al- 
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